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ABSTRACT 


It is known from the Apollo magnetometer experiments that 
significant electromagnetic induction takes place in the lunar in- 
terior. This induction is excited by fluctuations of the interplan- 
etary magnetic field and is detected by the induced fields on the 
surface of the Moon. This paper reviews these results breifly 
and discusses the formal properties of the theory. It is shown 
that the mathematical treatment parallels that for classical electro- 
magnetic scattering. Further the wavelength spectrum of the fluct- 
uations of the interplanetary magnetic field include scales con- 
sistent with the radius of the Moon. The consequence is that the 
Moon is excited in several modes. Quadrupole and possibly octu- 
pole magnetic multipoles are found in the data. The electric type 
radiation corresponding to transverse magnetic (TM) excitation 
appears suppressed and so far below the detection threshold of 
the magnetometers, but electric multipoles should be present 
as well. If so they would tend to be restricted to the outer layers 
of the Moon. 
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INTRODUCTION 

The purpose of this paper is to demonstrate the formal similarity 
between classical Mie scattering as known in optical phenomena and 
electromagnetic induction in the Moon caused by fluctuations of the in- 
terplanetary magnetic field. It is by now well known that the inter- 
planetary magnetic field contains a continuum of disturbances which 
are seen by the Moon as time dependent variations in the direction and 
strength of the magnetic field (Sonett, et al. , 1971). In turn these fluct- 
uations result in the flow of eddy currents in the lunar interior whose 
magnetic field is detected on the surface of the Moon by comparison of 
records taken there using the Apollo surface magnetometers and the free 
stream solar wind field simultaneously monitored by lunar satellite Ex- 
plorer 35 (Sonett, et al. , 1972). From the concurrent Fourier spectra 
of the surface and solar wind magnetic fields the internal bulk electrical 
conductivity profile of the Moon can be determined using the classical 
equations of electromagnetic scattering. 

The actual theory of scattering as is known for optical wave lengths 
and for vacuum conditions must be appropriately modified for the boundary 
conditions which exist at the contact between the Moon and solar wind. How 
ever the general formalism remains intact and we shall show later that the 
wavelength regime found in the solar wind which is responsible for the in- 
duction in the Moon is identifyable with magnetic multipole radiation including 
at least quadrupole and octupole terms. Proof of the effect of the solar wind 
in modifying the boundary conditions comes from observation that the lunar 
induction signal rises to values of 3-4 at the peak frequency, some 6-8 
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times that for vacuum induction. The increase is due to the dynamic 

pressure of the solar wind which forces the induced magnetic field 

lines back under the surface of the Moon. Skin depth arguments show 

that the confinement of the field varies with depth in the Moon; the 

confinement is frequency dependent and the shell within which the field 

is confined shrinks as the frequency is increased. 

The topology of the idealized interplanetary magnetic field lines is 

given by Parker where the equation for the lines is that of an Archimedean 

spiral (Parker, 1963). This follows from a combination of the solar 

rotation period, the high electrical conductivity of the solar gas, and 

the assumption of constant outward velocity for the solar wind. The 

latter is substantially correct for most of the flow field because of the 

strongly supersonic flow. Components of the Parker field (in the ecliptic 

2 

plane) are B = B (r / r) and B, = (cor/v )B where B is the component 
r r ro o t s r r r 

of the field radially outward from the Sun at r, B rQ is the radial component 

at the surface of the Sun, r is the solar radius, r the radius vector to 

o 

the position of observation, B^ the tangential component, w the angular 
velocity of the Sun, and v g the bulk flow speed of the solar wind. The 
angle of the field measured between the outward pointing radius vector 
and the tangent to the field line at the position, r , is = tan * (cor/v g ). 
Near Earth during quiet times, the ideal angle based upon v g = 400 km/ sec, 
is 45 deg. Thus the k vector field at the position of the Moon during 
quiet times is pointed at about 45 deg. away from the radius from the Sun. * 
During times of solar disturbance the field can take on any direction, and 
the spiral geometry is destroyed (Colburn and Sonett, 1966). Observation 
of induction by the solar wind field is limited to times when the Moon is in 

Footnote 


*Based upon Alfven waves propagating along the field lines. 
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the solar wind free of the magnetic tail of the Earth, though observation 
can be continued during those times when the Moon dwells in the post 
shocked plasma of the transition region behind the bow shock wave. 

The reason for the latter is that at lunar distance the post shocked plasma 
has returned nearly to the parameters found in the solar wind, and there - 
for the two are not easily distinguished. 

It is our aim in the following sections to show the formal similarity 
between the problem of lunar electromagnetic induction and Mie scatter- 
ing. For this demonstration we shall draw heavily upon earlier work. 

The key to the similarity lies in the applicable wave equations expressed 
in the next section which are seen to be just those of classical vacuum 
scattering of light from spherical objects. In the lunar case complica- 
tions arise from the radial stratification of the conductivity function for 
the Moon, modification of the boundary conditions because of the pressure 
of the solar wind which results in the formation of a surface current layer 
(Blank and Sill, 1969; Schubert and Schwartz, 1969), and the existence, 
in principle, of steady state excitation corresponding to ohmic dissipa- 
tion for the TM (transverse magnetic) mode which has no analog in vacuum 
scattering (Sonett and Colburn, 1967). Further for the Moon there is van- 
ishingly small radiation field at least for frequencies considered here; the 
induced field is primarily due to magnetic multipole excitation and comprises 
only a near field. Electric multipoles are by no means ruled out but have 
so far not been detected. The dominance of magnetic excitation arises 
naturally from the large fulk electrical conductivity. 

The basic equations written out in the next section have general applic- 
ability even for the more exact asymmetric induction problem which takes 
into account the existence of the lunar plasma cavity, but this paper is re- 
stricted to spherically symmetric induction. The two approaches are 
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not basically different, but asymmetric induction has been solved so far 
only for the low frequency limit (Schubert, et al. , 1972). Modification 
of the mode spectrum by the introduction of higher frequencies is presently 
being worked out and the complete solution should be available at a later 
date. 


FIELD EQUATIONS 

The TE (transverse electric) and TM (transverse magnetic) wave 
equations are respectively 

sj 2 J\ m + k lfL m = 0 d) 

and 

V 2 JI e -tL- ^ (/Xi e ) + = o (2 ) 

where equation 1 represents the transverse magnetic (TM) mode, equa - 

2 2 

tion 2 the transverse electric (TE) mode, k = E +£ f f*- co , 

kj = i u £ - and A* and E are the magnetic permeability and electrical 
permittivity respectively. The field vectors E and H are determined 
from the potentials and • The potentials and 

are given by 

CX: 

S\ m = /"V S H 0 sm? gA'J G™ <r) (cos 0 ) (3) 

1 


and 

n 


<iC' 


e j t a 

= H — cos 
o r 



fyA® G\ (r) (cos e ) 


where v is the speed of the solar wind, a the lunar radius , r, G , 
s 

spherical polar coordinates, (cos 6 ) the associated Legendre 


( 4 ) 



5 


polynomials and 


.fi 


a _ i (2^ + 1) 
V.1 ' T(F+T) 


The radial functions are solutions of the differential equations 


d 2 G 


m 


A + U 2 - -M+n. 


dr 


Cj 1 (r) - 0 


(5) 


and 


d 2 G| 

dr 2 


1 d(k 2 ) d G 1 


2 dr dr 


+ . A k 


§ 

2 . <io_t.iL ^ (r) . 0 (6) 


A more complete theory of the induction has recently been given by 
Schubert, et al. , (1972) which includes the asymmetry introduced by the 
flow field of the solar wind. However the fundamental properties of the 
induction can in most part be seen from the behavior of the spherically 
syymetric equations; departures are important but the latter theory is 
still under development for the time dependent case including the higher 
order time multipoles. 

Induction in the Moon is expressed most conveniently by the formalism 
of Schubert and Schwartz (1969) where 


H 
— s 

H 


= a 


& + <*■„/* VW “-’W , 


!>W 

V Kn fc 2 U.< k c R c 



cr B Pi<k c R c ) 
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( 7 ) 


for the magnetic field on the surface of the Moon. Here the first term 

represents the unit forcing field, the second the induction due to the TM 

mode and the last term that due to the TE mode. H q is the forcing field, 

a a unit vector in the direction of H , "X and ^ unit vectors in the 

— y i o— — 

directions along the axis of toroidal symmetry and poloidal symmetry 
respectively, f)” and (J”” are the shell and core conductivities for the 
model here of a two layer Moon, R and R the core and lunar radius 
respectively, JA- the magnetic permitivity, v g the solar wind speed, 
p^ and are spherical Bessel functions in the notation of Stratton 
(Stratton, 1941 ) and is an induction number given by 


B, (k R ) = 1 
1 ' c c 



3 cot(k R \ 
> c c' 


k R 
c c 


Lastly, 3 and cZ, are respectively the polar angles between the toroid- 
al or poloidal axes to the position of observation on the surface of the 
Moon and k^ is the wave number in the core Equation (7) represents 
the lunar response for the lowest order induction, but higher orders 
are obtainable from the general formulae of Schubert and Schwartz (1969). 
For illustrative purposes the dipole order suffices, until later when higher 
order TE interaction is specifically addressed. 

In equation (7) the second (TM) term is linearly dependent upon the 
speed of the solar wind, v , which determines the interplanetary electric 

field E = v x B which drives the TM mode. However, experimental 
— m -s — 
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findings show that this mode is suppressed and the lunar response is 
dominated by the TE mode as expected for a partially conducting 
material. We note that B, ( 1 for all k R and since R C R 
the denominator can never vanish; however, large values for the 

i 

ratio can be attained as (R /R m ) — f- 1 confirming that a strong 

current layer exists on the sunward side of the Moon consistent with the 
dynamic pressure of the solar wind forcing into the Moon the induced 
magnetic field and yielding amplifications well in excess of those 
found for vacuum induction. 

Although the customary treatment of scattering pays most attention 
to the far field, in the present case little radiation field is found. The 
near field dominates the response, and most attention is concentrated 
upon the interior of the scatterer, the Moon. We shall return to the 
properties of equation (7) in a later section. 

THE GEOMETRY OF LUNAR INDUCTION AND THE 
RESPONSE (TRANSFER) FUNCTION 

Viewed as a scattering problem the basic geometry of the excita- 
tion of lunar response by the solar wind electromagnetic field is shown 

$ 

in figure (1). E , B , and k are taken mutually orthogonally as is the 
case for Alfven waves propagating along the mean interplanetary mag- 
netic field which is taken along k (not shown). For the case of Apollo 12 
(from which the data shown in this paper are taken) the magnetometer lies 
nearly upon the lunar equator; the diagramatic representation of _k_ in 
figure (1) is indicative of the mean geometry but time dependent devia- 
tions take place which are significant and must be accounted for. Such 



details would not materially modify the present discussi on. The solar 
wind vector is shown as V making an angle % with k and V (phase 

s — p 

velocity). The position of the magnetometer (LSM) is also shown and 
© is the scattering angle measured from the direction of k to the 
outward pointing radius vector from the center of the Moon to LSM. 

The response of the Moon for modal order higher than dipole is asym- 
metric and qualitatively indicated by the two functions A & and A^. 

Data taken in the experiment discussed here are from the lunar 
satellite Explorer 35 which monitors the free stream solar wind mag- 
netic field free of perturbing influences of the Moon, and the LSM 
which determined the surface field on the Moon. 

Data taken in orbit are free of influence of the Moon because of the 
supersonic streaming of the solar wind which prevents upstream propa 
gation of waves from the Moon. 

We define an operational transfer function by 

h 2 .(f) + h u (f) = A(f). • h u (f) 

where h^^ (f ) > hj^(f) are the Fourier transformed ith components of 
the forcing field and response field, f the frequency, and A^(f) the 
ith component of the transfer function, all defined in a rectangular 
coordinate system where x is vertical locally at the Lunar Surface 
Magnetometer (LSM) site, y locally eastward, and z northward. 
LSM measures the sum of the forcing field and the induction while 
Explorer 35 in lunar orbit measures only the former. 

To give the reader some idea of the properties of the data obtained 
figure (2) shows the transfer function as defined above. These data 
have been corrected for plasma noise on the lunar surface, and are 
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representative of the response in the frequency interval shown. As 
we shall see shortly the rollover at high frequency and of the scale 
is indicative of excitation of higher order magnetic multipoles. The 
transfer function attains values higher than 3 verifying that containment 
of field lines is taking place. The model Moons shown in the insert 
are best fit monotonic models. They are used to make the forward 
calculated curves shown superimposed upon the data, and correspond 
to different values of phase velocity which is explained subsequently 
(See Sonett, et al. , 1972 for details). 

PARTIAL WAVES 

Equations (1) and (2) show that the partial wave excitation is 
practically restricted to the lowest order mode when the wavelength 
of the exciting radiation field is large compared to the scale size of 
the Moon, i. e. , X R m - This corresponds in the optical case to 
the Rayleigh limit. Actually for the Moon the problem is considerably 
more complex as higher order modes have been inferred in the itera- 
tion of the conductivity profile. 

Assuming the idealized case of Alfven waves propagating along the 
direction of the mean solar wind magnetic field, B , since B is 
spiralled because of the high conductivity of the plasma and the rota- 
tion of the Sun, the angle with respect to the radius vector from the 
Sun is about 45 degrees at 1 AU. Therefore the wave normals for the 
Alfven fields are slanted at about this angle. Since the angle of the mean 
field tan ^ ( r / v s ) changes depending upon the regional value of v . 

Other effects can also change the direction of B q over more or less short 
periods of time. Furthermore the wave field is not necessarily purely 
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Alfevenic all the time, and the direction of the wave normals, K , can vary. 

In general, however, a reasonable representation is obtained by the assumption 
of a mean field direction of about 45 deg. 

Figures 3 and 4 show the partial fields for two cases calculated 
using the best three layer Moon conductivd:y profile. Figure 3 is for a 
scattering angle © = 150 deg. and Figure 4 for © = 120 deg. Both use a 

phase velocity, v^ = 200 km/sec. somewhat lower than the usually value, but 
not seriously in error. The total induced field is given by , while the 
partial contributions are shown separated into real and quadrature parts. 

Both figures show significant increases in the higher order contributions as 
2 IT r / ^ increases, and the general conclusion that the Mie scattering 
range is attained is supported. 

Although not a central issue in this paper we show in Figure 5 all the model 
conductivity functions determined by iterative techniques ( Sonett , et al. , 1972). 
These indicate generally the magnitude of the internal conductivity and justify 
the conclusion that displacement currents are ignorable generally in the Moon 
except perhaps near the surface. The various models give a concise statement 
of the limitation of uniqueness permitting various combinations of monotonic 
and current layers.' All are electrically equivalent insofar as the amplitude 
of the induced field is considered, but differences in phase are likely. For 
the demonstrative purposes of this paper we use the best fit three layer (3L) 
model. 

Variations in v^ and in scattering angle lead to different scattering 
amplitudes when the latter is calculated using the 3L model. Several cases 
are shown in Figure 6 for comparative purposes. The use of a small value of 
v shifts the wave-length spectrum to smaller values and therefore increase 
the contributions from the higher order modes as is seen by the steeper 
rollover. Cases calculated are for a scattering angle of B = 150 degrees 
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which shows the rollover strongly. That the rollover is a strong function 

of scattering angle is seen by inspection of the upper part of the figure which 

uses a fixed value for v = 200 km/ sec. This clearly is an effect of 

P 

mutual interference between the different multipoles of the scattered field. 


WAVELENGTHS 

So far we have shown that the wave equations for lunar scattering are 
identical to those for the generalized case of induction in a spherical object; 
it is necessary to show that the wavelength regime fits that required for 
Mie scattering. The wavelength of disturbances in the solar wind is inferred 

V' 

from the Fourier transformed ,e£ time series record for each component 
of the magnetic field. This forcing field is measured on Explorer 35 lunar 
orbiter. Without knowledge of the solar wind speed it is impossible to 
infer the wavelength spectrum. Furthermore the wavelength is that taken 
normal to the surfaces of constant phase which generally are not oriented 
with k vector along the solar wind direction but along the spiral direction. 

We assume a mean velocity for the solar wind, v , = 400 km/ sec. , 
which is the approximate observable mean value. Next the direction of 
the mean interplanetary field is determined using minimum variance 
techniques. This is based upon determining the total spectral variance 
for each of the three components of the field in an arbitrary coordinate 
system and then rotating until the direction of minimum variance is found. 

The mean field direction is also found using an averaging technique. Analysis 
shows the direction of minimum variance and mean field to be colinear; 
discussion is restricted to these cases and they form the basis for the 
comments made here. 

Assumption of a reasonable solar wind speed and knowledge of the 
direction of the interplanetary magnetic and k fields permits a simple 
geometrical construction to be made to find the phase velocity and from 
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this the wavelength spectrum. Figure 7 shows a diagram of the geometry. 
The figure is a view down upon the plane of the ecliptic. The line OA is 
the radius vector from the Sun along which v g is measured. (It is assumed 
that the solar wind is free of aberration, unimportant in the present context. ) 
The line WW is the intersection of the ecliptic with the wave front; W^W is 
the wave front transported through the distance pp' by solar wind convection. 
(The propagation speed is small compared to this velocity. ). The wave- 
length is defined by QP' so that a wave front is moved normal to itself in 


the time it takes for the wave front WW to convect through pp'. Therefore 


the phase velocity, v = v cos 6 , and the wavelength is given by 

p s 

\ . where f is measured frequency. The 

wavelength regime for the frequencies encountered in the Apollo magnetometer 
experiment range from 7, 7X>\? ^ ^ ^ ^ measured in units of the 

scattering parameter ^ / A where R m is the lunar radius. Thus 

it is seen that the scattering range involved encompasses the Rayleigh 
limit at low frequency but extends well into the Mie regime at the high 
frequency end of the spectrum. The important consequence of this is the 
addition of higher order modes to the induction spectrum; the dominant 
magnetic multipole radiation of the excited field includes quadrupole, 
octupole and perhaps higher orders depending upon the forcing field con- 
figuration at the time of observation. Waves in the solar wind which are thought 
to be responsible for the bulk of the induction are consistent with Alfven 
waves propagating along the lines of force of the mean field. Such waves 
have been demonstrated by Belcher and Davis to be outward propagating and 
therefore "blue" shifted by the solar wind (Belcher and Davis, 1971). Their 
real frequencies are therefore likely still lower than those found by direct 
measurement in the frame comoving with the Moon. Such waves are known 
to be non-dispersive and have propagation speeds given by V a = B/(41^ ^ ) 
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where B is the field intensity, and ^ the mass density of plasma. For 

typical solar wind conditions V a 50 km/sec. , and therefore the prop- 
agation speed can be ignored in a first order analysis when compared to 
the convection speed, V gj of the solar wind. 

SPECIAL PROPERTIES OF THE MOON, POLARIZATION, ETC. 

This paper is restricted primarily to discussion of the electromagnetic 
excitation of the Moon by the Solar wind and formal similarity to Mie. 
scattering. We have discussed the excitation of higher order modes and 
shown this to be consistent with magnetic type radiation, but a radiation 
field is likely to be vanishingly small. Thus the far field in the sense of 
radiation theory is important, and the principle response of the Moon which 
is detected by magnetometers consists of the near field which is non- radiative. 
The conditions are analogous to an antenna at very low frequency where radi- 
ation is trivially small and only the close-in non- radiating electromagnetic 
field is seen. 

The Moon is a somewhat special object. Its high internal conductivity 
is consistent with the suppression of electric type excitation. As is well 
known such media are anamolously dispersive since the phase velocity 
increases with frequency. Displacement currents are ignored since 

where J is the bulk electrical conductivity, CjO the 
angular frequency, and the inductive capacity (all units and definitions 
in the S. I. units). Actually near the surface of the Moon such currents are 
expected but cannot be large and probably are below detectability. 

The TM mode can be shown by other arguments to probably be un- 
important, at least so far as analysis of existing data so far undertaken 
can be used as a guide. As pointed out earlier in this paper, the TM mode 
has anomalous behavior because of the presence of the conducting solar wind 
which permits ohmic currents to flow in the steady state. This has no analogue 
in vacumn scattering. It can be shown that the ohmic currents must decrease 
monotonically as frequency is increased because skin depth effects come into 
play and force the lunar interior current streamlines nearer to the surface. 
Therefore the total resistance increases with an attendent decrease in the 
integrated current. 

However the above argument cannot be correct for the reactive currents 
which flow in the TM mode. Using a simple two layer model for illustration, 
the capacity of the Moon can be shown to be the same as for a parallel plate 
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c apacitor of the same surface area. However the effective plate separation 
depends upon the skin depth assumed which depends upon the frequency. The 
plate separation will decrease as the frequency is increased, so that the 
capacitive current will increase with frequency. Therefore a peak in the TM 
response can take place at some frequency which depends upon various factors 
such as the conductivity profile in the lunar interior. 

Calculation suggests that these effects must be small; the absence of 
a detectible bow shock wave about the Moon shows that the steady state 
interaction is virtually suppressed. For these reasons present consideration 
of the Moon usually ignores the TM mode. 

Polarization of signals from the Moon has not so far been detected, but 
theory indicates that the induced field should be elliptically polarized. The 
reason for the difficulty is that the incoming Alfven wave field carried in the 
solar wind appears, when tested, to have random polarization. Until 
sufficiently "clean" signals, i. e. swaths of data with linearly polarized waves 
having fixed polarization are found in the overall data set, it is difficult if 
not impossible to make an adequate polarization test of the induction. Diffi- 
culties also exist wi th the coherence between Explorer 35 and the lunar surface 
because a decoherence by the intervening solar wind. Nevertheless theoretical 
calculations show that the expected lunar echo at the high frequency end of 
the data should be elliptically polarized. 


I 
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CAPTIONS 

1. Geometry of the scattering process. The plane containing E and B is 
the plane of the disturbance wave in the solar wind. The propagation 

is shown along k which is also the direction of the interplanetary magnetic 
field for Alfven waves. V g is the solar wind bulk velocity and v^ the phase 
velocity including the convection of the disturbance which dominates the 
actual propagation of waves in the solar wind. The position of the surface 
is indicated by LSM and the scattering angle, 0 , is defined as the angle 

between k and the position vector of LSM referred to the center of the Moon. 
Aq and A ^ are the response functions in the east-west and north- south 
directions; they are different for higher order modes. 

2. Transfer function amplitude for lunar response corrected for plasma 
generated noise at the lunar surface. The insert shows best fit model 
conductivity profiles obtained by inversion of the transfer function together 
with an assumed value for phase velocity, v^. Forward calculations of 
the lunar response are shown as the lines superimposed upon the data. 

3. Partial wave contributions to the scattered (near) wave field showing the sum 

( ), and the real and imaginary parts of the partial waves for^, ^1^5 

assuming v^ = 200 km/ sec, and scattering angle ® - 150 deg. 

4. Similar to Fig. 3 except for scattering angle & = 120 deg. 

5 o Lunar internal conductivity profiles showing two layer (2L), three layer 
(3L), four layer (4L), core plus layer (CCL), current layer (CL), and 
double current layer (DCL) cases. All are approximately equivalent 
electrically, but significant differences do exist in the residual errors. 

Also no account has been taken of phase up to this time. For details 
see Sonett et al, (1972). 
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Forward calculation of transfer functions based upon the best fit three layer 

model of Fig. 5. The effects of varying v^ while holding & fixed are shown 

in the lower part of the figure. 'All cases shown some rollover. The upper 

part of the figure demonstrates the effect of changing scattering angle while 

holding v fixed. 

P 

Geometry of plane waves propagating and convecting with the solar wind. 

The wave front WW moves one wavelength normal to itself in the time 
taken for the distance PP' tp be covered by solar wind convection. Thus 
wavelengths must be corrected downwards by the factor cos & . 
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